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The ISO-LWS map of the Serpens cloud core: 
I. The SEDs of the IR/SMM sources* 
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Abstract. Iso-Lws mapping observations of the Serpens 
molecular cloud core are presented. The spectral range is 
50 — 200/im and the map size is 8' x 8'. These observations 
suffer from severe source confusion at Fir wavelengths 
and we employ a Maximum Likelihood Method for the 
spectro-spatial deconvolution. The strong and fairly iso- 
lated source SMM 1/FIRS 1 presented a test case, whose 
modelled spectral energy distribution (SED), within ob- 
servational errors, is identical to the observed one. The 
model results for the other infrared and submillimetre 
sources are therefore likely to represent their correct SEDs. 
Simulations demonstrating the reliability and potential of 
the developed method support this view. 

It is found that some sources do not exhibit significant 
Fir emission and others are most likely not pointlike at 
long wavelengths. In contrast, the SEDs of a number of 
SMMs are well fit by modified single-temperature black- 
bodies over the entire accessible spectral range. For the 
majority of sources the peak of the SEDs is found within 
the spectral range of the Lws and derived temperatures 
are generally higher (> 30 K) than have been found by ear- 
lier deconvolution attempts using Iras data. SMM sizes 
are found to be only a few arcsec in diameter. In addi- 
tion, the SMMs are generally optically thick even at Lws 
wavelengths, i.e. estimated A(t = 1) are in the range 160- 
270 ^m. 

The Rayleigh- Jeans tails are less steep than expected 
for optically thin dust emission. This indicates that the 
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SMMs are optically thick out to longer wavelengths 
than previously assumed, an assertion confirmed by self- 
consistent radiative transfer calculations. Models were cal- 
culated for five sources, for which sufficient data were 
available, viz. SMMl, 2, 3, 4 and 9. These models are 
optically thick out to millimetre wavelengths (wavelength 
of unit optical depth 900 to 1 400 /ini) . Envelope masses 
for these SMMs are in the range 2-6 Mq, which is of 
course considerably more massive than estimates based 
on the optically thin assumption. The luminosities are in 
the range 10-70 Lq, suggesting the formation of low-mass 
to intermediate mass stars, so that the existence of such 
massive envelopes argues for extreme youth of the SMMs 
in the Serpens cloud core. 

Finally, we present, for the first time, the full infrared 
SEDs for the outburst source DEOS, both at high and low 
intensity states. 

Key words: ISM: individual objects: Serpens cloud core 
- clouds - general - dust, extinction - Stars: formation 



1. Introduction 

The Serpens molecular cloud is a magnificent laboratory 
for the study of multiple low mass star formation. The 
cloud core has a high density of young stellar objects 
(YSOs), including several submillimetre sources: for a dis- 
tance of 260 pc, Kaas (1999) estimated a surface density of 
400 - 800 YSOs per square parsec, which is several times 
higher than found in the star forming western parts of the 
pOph clouds (~ 100 pc~^, Bontemps et al. 2000 and refer- 
ences therein). The extremely high visual extinction makes 
the Serpens cloud core nearly exclusively accessible to in- 
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Table 1. ISO observations of the Serpens cloud core 
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Fig. 1. Outline of the 8' x 8' map of the Lws observa- 
tions, shown together with the 450 /xm continuum map of 
Davis et al. (1999). Conventional source names are indi- 
cated. Solid circles about grid-points and broken circles 
about other map-points (see Table |l|) outline the contours 
of the Hpbw of the Lws, taken as 70". The small sym- 
bol south of SMM9 identifies the pointing of the Sws, the 
aperture of which is rectangular 14" x 20" to 17" x 40" 
depending on the wavelength. The Lws-map is centred on 
the far-infrared/sub-millimetre source FIRS 1/SMM 1 and 
the offsets of 100" in both Right Ascension and Declina- 
tion (plus signs) correspond to the pre-flight estimates of 
the Lws-beam 
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Mode''^ 


Offsets 
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Date 


Name 


Investigator 




R.A. 


Dec. 




LWS 


Map 


R. Liseau 


LOl 








21-Oct-96 


Flow 


P. Saraceno 


LOl 
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23-Oct-97 


Field 


P. Saraceno 


LOl 


-100 


50 


08-Mar-97 


Field 


P. Saraceno 


LOl 


100 


50 


08-Mar-97 


SMM2 


M. Spaans 


L02 


151 


-148 


14-Apr-96 


SMM3 


M. Griffin 


LOl 


140 


-78 


21-Oct-96 


SMM4 


M. Griffin 


LOl 


101 


-124 


21-Oct-96 


S68N 


M. Spaans 


L02 


-32 


95 


14-Apr-96 


SWS 


DEOS 
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T. Prusti 
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61 
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61 
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T. Prusti 
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S68N 


P. Andre 
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G. Olofsson 
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G. Olofsson 
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61 
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Notes to the table: 

^ Instrumental modes: LOl and SOI refer to the full wavelength 
range of the Lws and Sws, respectively. L02 stands for Lws 
line spectrum scans. CVF denotes spectrum imaging with the 
Continuous Variable Filter of ISOCAM. LW2 and LW3 are two 
broad band IsoCAM filters centered at 6.7 and 14.3 /im respec- 
tively. 

* Offsets are relative to the centre coordinates of the Lws map, 
viz. R.A. = 18''29'"50'29 and Dec. = 1° 15' 18''6, J 2000. 



frared and sub-/millimetre observations (e.g., Casali et al. 
1993, Hurt & Barsony 1996, Davis et al. 1999). The sub- 
millimetre sources, named SMM by Casali et al. (1993), 
are distributed in a southeast-northwest direction with a 
concentration of sources in two clusters. This is also true 
for the presumably more evolved objects (Class I and II, 
see: Kaas 1999). The region also includes a great number 
of Herbig-Haro objects and molecular flows, where at least 
some can be connected with SMMs (e.g., Ziener & Eisloffel 
1999, McMulHn et al. 1994, White et al. 1995, Curiel et 
al. 1996). Outflows are seen to accompany star formation 
at all observed evolutionary stages. 

Submillimetre sources are often found in dense molec- 
ular cloud cores and are believed to identify birth places 
of stars. For the Serpens sources, Hurt & Barsony (1996) 
used image sharpening techniques to analyse low resolu- 
tion, broadband Iras data. Cores in low mass star forming 
regions are generally thought to be cold, with tempera- 
tures of about 20 K or below. Emission from such objects 
will peak at wavelengths typically longer than 100 /im, 
which makes the Long Wavelength Spectrometer (Lws) 



on board the Infrared Space Observatory (Iso) an instru- 
ment more suited than Iras to study these sources. In ad- 
dition, the Lws provides two orders of magnitude higher 
spectral resolution over its spectral range of 50 to 200 ^m. 

The Lws observations presented in this paper aim to 
address the nature and evolutionary status of the SMMs 
in the Serpens cloud core, in particular their physical state 
(T, p, M, L etc.). The observables to be discussed are 
the spectral energy distributions of these objects, obtained 
over a broad spectral range. This paper utilises mapping 
observations of the Serpens cloud core with the Lws and 
complementary observations, including pointed Lws data 
and observations with other instruments aboard ISO. In 
Sect. 2, an account of these observations (8' x 8' map and 
Iso archive retrievals) and their reduction is given, and 
the results are presented in Sect. 3. The Iso-Lws map is 
highly undersampled and in Sect. 4, we introduce a Max- 
imum Likelihood Method for the extraction of the SEDs 
of the spatially confused submm sources. The discussion 
focusses on the comparison of conventionally applied anal- 
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ysis methods (single temperature, optically thin approxi- 
mation) with that based on the detailed and self-consistent 
modelling of the radative transfer. Finally, in Sect. 5, the 
main conclusions of this work are briefly summarised. 

2. Observations and data reductions 

2.1. The ISO observations 

A grid of 25 positions of FiR-spectra in the 
Serpens cloud core was obtained with the Long- 
Wavelength Spectrometer (Lws; 43 - 197 /xm, 
R\ = 140 — 330) on board the Infrared Space Ob- 
servatory (ISO) on October 21, 1996. An account of the 
ISO-project is given by Kessler et al. (1996). The Lws 
is described by Clegg et al. (1996) and Swinyard et al. 
(1996). 

The formal map centre, viz. a — 18'^29™50^29 and 6 = 
1° 15' 18'.'6, epoch J 2000, coincides to within 10" with the 
position of the sub- millimetre source SMMl. The point- 
ing accuracy in the map is determined as 1" (rms). The 
spacings between positions in the map, oriented along the 
coordinate axes, are 100" in both Right Ascension and 
Declination. These offsets correspond to the pre-flight esti- 
mates of the circular Lws-beam, Hpbw. The actual beam 
widths are however smaller, 70" to 80", leading to appre- 
ciable spatial undersampling. This is also evident in Fig.|l|, 
where the Lws-beams and pointings are represented by 
circles. The size of the 5x5 Lws-map is thus 8' x 8', cor- 
responding to (0.7 X 0.7 = 0.5) I?3jQ pc^, where D^io de- 
notes the distance to the Serpens cloud core in units of the 
adopted value of 310 pc (deLara et al. 1991). 

At each map-point the grating of the Lws was scanned 
6 times in fast mode, oversampling the spectral resolution 
by a factor of 4. Each position was observed for nearly 
15min. The centre position was re-observed half a year 
later on April 15, 1997, for a considerably longer integra- 
tion time (24 spectral scans; Larsson et al. in prep.). The 
internal agreement between the deep integration and the 
map spectrum is excellent, giving confidence that the rest 
of the map data is also of good quality. 

Our map spectra were complemented with Lws spec- 
tra at 7 positions inside the mapped area, which were re- 
trieved from the Iso-archive (Table Q). In addition, ID- 
spectral observations with the Short- Wavelength Spec- 
trometer (deGraauw et al. 1996: Sws; Fov=14" x 20" 
to 17" X 40", 2.4 - 45 ^m, Rx ~ 10^ - 10^) and imag- 
ing spectrophotometry with the Continuous Variable Fil- 
ter (Cvf; 2.5 - 15.5 ^m, pixel-Fov = 6", Rx > 35) of Iso- 
CAM (Cesarsky et al. 1996) were also analysed. Further, 
broadband ISOCAM images at 6.7 and 14.3 /im and with 
pixel-Fov both 1"5 and 3" were part of the observational 
material. 

In summary. Table |l| provides an overview of the ob- 
servational material employed in this paper. Fig. ^ shows 
the Sws and Lws pointing positions superposed onto the 
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Fig. 2. 50 - 200 /im spectra towards the grid points of the 
Lws-map. The spatial offsets are given along the upper 
and right-hand axes. The wavelength scale is found along 
the lower axes and that of the flux density along the left- 
hand axes, except for (0", 0") and (-100", -100"), for 
which the scale is given inside the sub-frames. The Lws 
data have been rebinned to the spectral resolution Rx — 
20. The error bars reflect the differences of the displayed 
median data (solid line) and the original data, prior to 
the detector stitching. The FiR-background spectrum is 
defined as the straight average of the spectra inside the 
boxes with the thick lines 

sub-millimetre continuum map (450 /im) of Davis et al. 
(1999), where 9 discrete sub-millimetre sources (SMMl - 
SMMIO, except SMM7) are identified. 

2.2. The reduction of the ISO data 

The Lws observations were pipeline-processed (0LP7.1) 
and were subsequently reduced with the interactive analy- 
sis package LIA. Post-pipeline processing of the Sws and 
Lws data was done with the package ISAP and of the 
ISOCAM data with the corresponding CIA programs. 

At each position, the individual Lws scans for each of 
the 10 detectors were examined, 'deglitched' and averaged. 
Corrections were applied to the 'fringed' spectra in the 
map. The fringing indicates that the emission is extended 
and/or that point sources were not on the optical axis of 
the Lws, i.e. the radial distance from the optical axis was 
larger than about 25". 

For the Lws data we believe that the absolute flux cal- 
ibration is good to an accuracy of 30% (Swinyard et al. 
1996), whereas relative offsets between overlapping spec- 
tral regions of adjacent detectors were generally within 
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Fig. 3. Spatial distribution of the continuum emission in 
the Lws-map in wavebands centred on 50, 60, 80, 100, 
130 and 160 /xm, with spectral bandwidth of 10 /im. Con- 
tours are spaced by 210~^^W cm~^, the lowest contour 
is at 2 10"^'' W cm~^ and the displayed dynamic range is 
a factor of 10 
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Fig. 4. A pure blackbody fit (/3 = 0, see the text) to Lws 
data of SMM 1 results in the indicated values of the tem- 
perature and diameter of a circular source, the formal er- 
rors of which are given as the standard deviation about 
the best- fit values (x^-minimisation). The displayed Lws 
spectrum has been re-binned to a resolution of 20. The 
displayed error bars correspond to the relative flux uncer- 
tainty, when stitching together the 10 individual detectors 
of the Lws. At the bottom of the figure, the (subtracted) 
average background spectrum is shown for reference and 
compared to Iras data (filled circles) 

10% ('detector stitching' uncertainty, see also Fig.||). In- 
ternally, the Lws-accuracy is much higher than 30% for 
the Serpens cloud core, as is evidenced by observations of 
the same targets at different times and under different 



observing conditions (e.g., pointed vs. mapping observa- 
tions). 

The absolute accuracy could be less, however, in the 
very weak spectra towards the northern and western edges 
of the map, for which the flux uncertainties are dominated 
by the dark current correction. This was estimated based 
on measurements directly before and after, but not during, 
the mapping observations. Suspected dark current drifts 
during the observation of the first 5 map points, viz. from 
(-200", "200") to (-200", 200"), were however approxima- 
tively linear in time, whereas for the remaining spectra, 
the dark current seemed to have stabilised at a constant 
value. Support for our applied dark current corrections is 
provided by the fact that these resulted in fluxes which 
arc in very good agreement with IRAS-ISSA data at 60 
and 100 /im (formally within 7% and 4%, respectively, see: 
Fig. I). 

The reduction of the Sws data is similar to that of the 
Lws data, whereas the reduction of the CvF data involved 
additional corrections for transients and flat field. 

3. Results 

From Fig. || it can be seen that detectable Fir emission 
is present at all observed positions. The sensitivity of 
these map data at long Lws wavelenghts is of the order 
of 10^^^ W cm^^ jjjnT^ (about 3 Jy per beam at 100 /xm), 
whereas this noise level becomes worse by a factor of about 
two towards the shorter wavelengths. 

As pointed out in Sect. 2.2, most of the observed spec- 
tra were fringed. The situation is different at the map cen- 
tre, where the Lws spectrum shows hardly any fringes at 
all. This strongly suggests that the FiR-emission is domi- 
nated by a source which is pointlike to the Lws and which 
was reasonably well centred in the aperture during the 
observation. As is also evident from Fig. |[ where Lws 
measurements in 6 continuum bands between 45 /im and 
165 /im are presented, that this source dominates the emis- 
sion at all Fir wavelengths. In the following section, the 
physical characteristics of this source will be discussed. 

4. Discussion 

J^.l. The sub-millimetre source SMM 1 

4.1.1. The blackbody observed by the LWS 

The centre position of the map dominates the emission 
over the entire wavelength range of the Lws (Fig.|^). This 
map point coincides with the strong sub-millimetre source 
SMMl/FIRSl. At Fir and sub- mm wavelengths, this 
source is relativly isolated. Any problem of source con- 
fusion is therefore very much reduced in this case. 

As evident in Fig. ^ the Fir spectrum is well fit by 
a single blackbody of temperature 36 K. Such a relatively 
high temperature was also determined by, e.g., Nordh et 
al. (1982), Harvey et al. (1984), McMuUin et al. (1994) 
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Fig. 5. The SED (A Fa vs A) of SMM 1, between 2 /im and 6 cm. Open symbols identify data which have been obtained 
with apertures larger than 7", with the small open symbols indicating the according to Eq. (3) beam-matched data. 
The thin solid line shows a conventional (modified blackbody) fit to these data within 40 /im J; A 3 mm and, for 
reference, the 36 K-blackbody of Fig. ^ is shown by the dash-dot line. The dashed Hne depicts the SED obtained from 
a model, for which the radiative transfer through the dust envelope is treated self-consistently. The line thickness for 
the Lws and CvF data, respectively, are shown in the legends inside the figure. The accompanying upper limit symbol 
refers to the short-wave CvF data. The photometric data collected from the literature are referenced to as K = Kaas 
(1999), HWJ = Harvey et al. (1984), HB = Hurt & Barsony (1996), DMRDR = Davis et al. (1999), CED = Casah 
et al. (1993), MMWHB = McMuUin et al. (1994), HDSB = Hogcrheijde et al. (1999), CRMC = Curiel et al. (1993). 
The near-infrared data (A < 10 /im) are probably unrelated to SMM 1 (see: the text) 



and Davis et al. (1999). The diameter of the hypothet- 
ically circular source is ~7", i.e. (2 200 ± 250)1)310 AU. 
This latter value confirms the point-like nature of the far- 
infrared source, as was suspected in the previous section 
on the basis of its largely unfringed spectrum. It is also 
comparable to the source sizes seen in the interferomet- 
ric observations at millimetre wavelengths by Hogerheijde 
et al. (1999). In contrast, the temperature derived here 
is significantly higher than the 27 K obtained by Hurt & 
Barsony (1996) from high resolution iRAS data (HIRES). 



4.1.2. A modified blackbody fit to the SED of SMM 1 

In Fig.|^, the Lws spectrum is shown together with data 
collected from the literature. From this overall Spectral 
Energy Distribution (SED) it is clear that the Lws data 
are sufficient to determine the luminosity of the source. 
However, it is also apparent that at longer wavelengths 
the SED is steeper than the Rayleigh- Jeans fall-off of the 
blackbody (A Fx oc A~^), which could indicate that the 
emission is becoming optically thin somewhere outside the 
spectral range of the Lws. Provided that in this optically 
thin regime the dust opacity follows a power law, viz. k\ oc 
A"'', this index (3 can be determined from the observa- 
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tions, since then A Fx oc A~(^+'^) as A ^ cx) for any density 
and/or temperature structure of the envelope (spherical 
symmetry and central radiating source). A least squares 
fit to the data between 0.8 mm and 3.4 mm yields /Sobs — 
- (^^^|x^ + 3) = 1.2±0.2. At these relatively long wave- 
lengths, this statistical error is considerably larger than 
any deviation of the modified Rayleigh-Jeans law from 
the (modified) Planck curve, viz. FA(Planck)/_F'A(R.J) — 
(1.44 X lOVA^n, T) I [exp (1.44 x 10VA^„, T) - l] , which 
for T = 36K would result in the "correction" of /3obs by 
-0.087. 

Following standard procedures (e.g. Hildebrand 1983, 
Emerson 1988), we can then estimate the envelope mass. 
Using only the data between 40 /im and 3 mm, which were 
collected with sufficiently large telescope beams, i.e. larger 
than 7" and shown as open symbols in Fig.|5[ a modified 
blackbody can be fit to this section of the SED for a single 
temperature. Hence, the fiux density is given by 



^x^sBx (T) 



(1) 



with obvious notations. In this wavelength regime, the 
opacities of dust grains generally display power law de- 
pendencies (e.g., Ossenkopf & Henning 1994). Hence, let 
the continuum optical depth scale as 



/3 



(2) 



A first order correction for the finite beam sizes to the 
flux density (shown as small open symbols in Fig.|^) was 
included by 



Fx 



A,LWS / ^\,S + ^\,B 



A, B 



A,LWS 



Fx,B 



(3) 



where JIa, lws is the beam size of the Lws at different 
wavelengths, A (see Section 4.2.1 and Fig.^, flx,B is the 
beam size of the telescope in question, Qx, S is the source 
size, and Fa, b is the source flux reported in the literature. 
Any dependence of the source size on the wavelength is 
naturally accounted for by the radiative transfer modelling 
of Sect. 4.1.3. Here, it is assumed that the source size is 
the same as that of the blackbody fit and constant with 
wavelength. 

This modified blackbody remains optically thick far 
beyond 100 /im, up to A(t = 1) = 250 /im, the wavelength 
of unit optical depth. Longward of this wavelength the 
emission is optically thin and the total column density, 
A^(H2), can be estimated from 



TA = J Kx p{t) d£:^ Kxfi mn A^(H2) 
so that 



gas 



dust 
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Fig. 6. CvF view of the 5 to 15 /im SED towards SMM 1 
inside the Lws beam. The pixel size is 6", hence the shown 
map is 42" x 30". The map orientation is offset from the 
axes of the equatorial sysrem, rotated at a negative po- 
sition angle of a few degrees. Note the clear difference 
between the SED of the NiR object (distinctly 'blue') and 
that of SMM 1 (distinctly 'red') close to the map centre 
and separated by less than 20" 



Taking kx,, ~ Icm^g^^ at Ao — 1.3mm (Ossenkopf 
& Henning 1994) and ta^ = 0.05, fi — 2.4 for molecular 
gas and the gas-to-dust mass ratio Mgas/Mdust = 10^, the 
average column density is A^(H2) — 1.3 10^'* cm~^. For the 
assumed source radius of 1100 AU, the mass of the opti- 
cally thin envelope of SMM 1 is 2.2 Dfn, Mq. The average 
volume density n(H2) corresponds then to several times 
10^ cm""^, which is consistent with our choice of the theo- 
retical opacities for grains with thin ice mantles for dense 
protostellar regions. However, the observationally deter- 
mined value of P seems only marginally compatible with 
the opacity law, for which f] = 1.8. If one accepts that the 
opacities are generally correct per se, the observed low /3 
values (see also Table ||) could indicate that the optical 
depth through the dust has been under-estimated even at 
long wavelengths. 

The assumption of an isothermal, homogeneous sphere 
is of course very crude and physically not at all compelling. 
However, the Lws data do not contain any information 
which can constrain the structure of the source. Hoger- 
heijde et al. (1999), on the other hand, found evidence 
for structure in the visibility curves of their mm-wave 
data and adopted a radial power law distribution for the 
density, n{r) oc with p ~ —2, so that their model of 
SMM 1 is more centrally condensed. For a spherical, cen- 
trally heated optically thin dust envelope, a power law is 
also expected for the temperature profile, viz. T(r) oc r''. 
The exponent is in this case given hy q — —2/(4 + /3), 
where /3 refers as before to the wavelength dependence 
of the absorption cross section of the dust grains. Since 
/3obs ~ 1, the temperature profile could be expected to 
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Fig. 7. The radial run of the temperature in the (assumed) 
spherical envelope of SMM 1. The solid line shows the opti- 
cally thick dust envelope, obtained from the self-consistent 
radiative transfer model discussed in the text (LLM = 
this paper). For comparison, the optically thin model with 
T oc r^" '' by Hogerheijde et al. (1999), and identified as 
HDSB, is shown by the dashed line. Both models assume 
an density law and the same grain opacities. When 
compared to the optically thin assumption, the run of the 
temperature is markedly different: a considerably steeper 
gradient is needed in most regions of the envelope to drive 
out the centrally generated flux (^ 70 Lq) 



exhibit q —0.4, a value which was indeed also used by 
Hogerheijde et al. (1999). 

It has been suggested that SMM 1 is a so called Class 
source (Andre et al. 1993) and the correct assessment of 
its envelope mass (and central mass, of course) is of par- 
ticular importance for such type of object. In the following 
section, we shall attempt to improve on the reliability of 
the parameters for SMM 1 by treating the radiative trans- 
fer through the dust envelope in a self-consistent manner. 

Before doing that, we need to address the possible 
near- infrared excess over a Class SED, however (see 
Fig.||). It is clear that over the 7" source size, the deduced 
column density, 7V(Il2) = 10^'' cm~^, would imply an ac- 
companying visual extinction of Ay ^10'^ mag (Bohlin et 
al. 1978, Savage & Mathis 1979), or still nearly Ak 
100 mag at 2 /im (Rieke & Lebofsky 1985). Therefore, on 
arcsecond scales around SMM 1 one would hardly expect 
any detectable direct flux in the near-infrared. This asser- 
tion is confirmed by Isocam-Cvf observations (Table |l]) 
which show that the NiR flux originates from a source 
different from SMMl. In Fig.|, a map of the spectra on 
the 6" pixel scale are shown. Obviously, SMM 1 and the 
source west of it do not have the same SEDs in the 3 to 
15 iim region. Hence, we conclude that the SED of SMM 1 
contains less NiR flux than indicated in Fig.^ and, conse- 
quently, these observations are supporting the view that 
SMM 1 is indeed a Class object. 
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Fig. 8. Surfaces of unit optical depth of the model of 
SMM 1 in the far infrared and sub-millimetre (40 /im to 
about 1 mm) . Shown is the variation of the source size 
with wavelength, i.e. the radius, R{t\ = 1) in AU, and 
the corresponding dust temperature at this radial distance 
from the central heating source, T{t\ = 1) in K. The dash- 
dotted line, labelled R{F\I2), designates the locus of the 
radius at which half of the flux at a given wavelength is 
emitted 



4.1.3. Self-consistent radiative transfer models of SMMl 

The models by Hogerheijde et al. (1999) are based on the 
assumption of optically thin emission. However, it was re- 
alised in the previous section that optical depth effects 
might affect the SED of SMM 1 even at long wavelengths. 
In order to test this hypothesis we have run models of 
the transfer of radiation through the dust envelope of the 
source. A fuller description of these computations will be 
presented in a forthcoming paper (Larsson et al. in prep.). 
From the comparison with the observations, these self- 
consistent models provide the temperature distribution in 
the envelope. For instance. Fig. ^ depicts the true temper- 
ature profile of the Hogerheijde et al. model of SMM 1 (see 
below). For the general cloud background we estimate a 
temperature of 20 K, and that flux has been subtracted 
from the model of SMM 1. 

Complying with Hogerheijde et al. (1999), we adopt 
a spherically symmetric geometry of the dust envelope 
about a central heating source and a radial power law 
distribution for the density, with p = —2 (in order to 
reproduce the visibilites obtained by these authors), i.e. 
n{r) — 1.7 X 10^^ (r/12 AU)~^ cm~'^, which reproduces 
the density at 100 AU we inferred from that paper. As be- 
fore, the opacities were taken from Ossenkopf & Henning 
(1994), specifically the model identified as MRN with thin 
ice mantles, density 10^ cm~'^ and age 10^ yr (although 
one would not expect these thin ice mantles to survive in 
the central regions of the envelope, where temperatures 
are greately in excess of 100 K, we kept this particular 
choice of opacities for comparison with the other models; 
using other grain compositions did not change the general 
results) . 
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Table 2. Dust temperatures, sizes and opacity indices of IR/SMM envelopes in the Serpens cloud core 



source iName 




-tdust (Kj 






(arcsec) 






Pobs 




Note^' 


Db ' 


1 1 


HB 


Db' 


1 1 


HB 


T T H 


DMRDR 


CED 


SMMl/FIRSl 


36 


43 


27 


7± 1 


4 


12 


1.2 ±0.2 


0.9 


0.5 ±0.2 




SMM2 


33: 


36 


24 


3 ±2: 


1.5 


20 


0.2 ± 1.9 


1.1 


0.8 ±0.8 


1 


OiViiVi o 




32 


24 


zt Z 


o 

O 


20 


0.6 ±0.2 


0.8 


1-1-0/1 




SMM4 


30 


27 


20 


5±2 


4 


20 


0.8 ±0.2 


0.8 


0.5 ±0.3 




SMM5/EC53 












20 


0.9 ±0.9 


0.9 


1.6 ±0.5 


1,2,3 


SMM 6/SVS 20 












20 


0.5 ± 2.4 


0.8 


1.3 ±1.0 


1,2,3 


SMM8 


21: 






6 ±3: 














SMM 9/S 68N 


24 


27 


23 


9±4 


3.5 


20 




1.2 






SMMll/HBl 


31: 






4±2: 












1,3 


HE 2 


24: 






5 ±3: 














SVS2 


25: 






9 ±5: 












1,3 


References: hb — 


Hurt & 


; Barsony (1996), llm 


— this 


paper, dmrdr 


= Davis et al. (1999), CED = 


Casah et al. (1993) 





Notes to the table: 

^ 1 = extended emission, 2 = excluded from calculation, 3 — near infrared source (see the text). 

^ Temperature and diameter values are for bb = blackbody fit and rioo^jm ~ 1 refers to the radiative transfer model. 

^ Our /Jobs has been evaluated between 0.8 mm and 3.4 mm. 



Beam Plateau Amplitude Width 




Fig. 9. The beam of the Lws at the centre wavelengths of its ten detectors (SW 1 to SW 5 and LW 1 to LW 5, see Clegg 
et al. 1996). Azimuthal averages are represented by the filled dots and the full lines refer to the continuous functions 
fitted to these data. For the distance r from the beam centre (optical axis) in arcsec we have il.\{r) = 1 ii r < P 
and 0,\{r) = A exp[— i (^)^] if r > P, where the plateau P is shown in the second panel from the left, the amplitude 
A in the following panel and the width tq in the right hand panel. An explanatory sketch of this parameterisation is 
provided in the farniost left hand panel 



The model reproduces indeed the visibilites of SMM 1. 
Further, the observed SED is reasonably well fit at long 
wavelengths (Fig. ||), i.e. the observed /?- value is obviously 
reproduced. This is mainly because the envelope is opti- 
cally thick out to nearly 1 mm (At-=i = 900 /xm, see the Ta- 
bles || and ^, where also other model results are provided). 
As expected, the "blackbody parameters" of the model are 
not vastly different from the single temperature fit: e.g., 
the diameter and temperature of the spherical source are 



4" and 43 K, respectively, for unit optical depth at 100 /xm, 
viz. for Tioo^^m = 1- Values at other wavelengths are pre- 
sented graphically in Fig.||. Evidently, at the longer wave- 
lengths, one 'looks' progressively deeper into the envelope, 
towards hotter and denser regions. Further, at the longer 
wavelenghts more than half of the flux originates from the 
inner regions, where the dust destruction front is located. 
This is shown in the figure as R{F\/2), i.e. the radius at 
which half of the flux at a given wavelength is emitted. 
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Table 3. Luminosities, opacities and masses of SMM envelopes in the Serpens cloud core 



Source Name 




(L 


0) 






iV(H2)* 
(10^4 cm-2) 




enve 


ope (M 


3) 


Note 


LLM 


HB 


CED 


I 


II 


I 


II 


HB 


HDSB 


SMMl/FIRSl 


71 


46 


84 


250 


900 


1.3 


2.2 


6 


3 


5.2 




SMM 2 


10 


6 


0.8 - 32 


170 


970 


0.7 


0.3 


2.3 


0.6 




Extended 


SMM 3 


18 


8 


2-32 


170 


1000 


0.7 


0.6 


4 


0.9 


1.8 




SMM 4 


16 


9 


3-32 


270 


1400 


1.5 


1.2 


6 


3 


3.2 




SMM9/S 68N 


16 


6 




160 


1200 


0.6 


1.7 


5 


1.0 







References: llm = This paper, HB — Hurt & Barsony (1996), CED = Casali et al. (1993), I — Single temperature/optically thin 
approximation, ll — Radiative transfer model, hdsb — Hogerheijde et al. (1999). 
Notes to the table: 

^ All reduced to the same distance, viz. d = 310 pc. Our luminosities refer to 40 - 3 000 fim. 

^ Opacities from Ossenkopf & Henning (1994: MRN with thin ice mantles, 10*^ cm~^ and 10^ yr, for which (3 = 1.8). 
* Refers to l (single temperature, optically thin approximation). 



We shall not dwell further on the details of these spe- 
cific radiative transfer calculations, for the following rea- 
son: although the model fit is quite acceptable in the FIR 
and submm/mm spectral region, it clearly underproduces 
the flux at wavelengths shortward of the peak of the SED 
(Fig.^. The remedy of this shortcoming requires the re- 
laxation of the assumption regarding the source geometry. 
Because of disk formation, spherical geometry is generally 
not considered a good approximation for YSOs, including 
protostellar objects. It is possible, in fact, to satisfactorily 
fit the entire SED using a radiative transfer model in a 2d 
axially symmetric geometry (disk plus envelope). The re- 
sults of such computations will be presented in a follow-up 
paper (Larsson et al. in prep.). 



4-.2. The other SMMs of the Serpens cloud core 

In Fig.|^, a secondary emission peak lies on the extension 
towards the south-east. This morphology is very similar to 
that previously seen at Fir wavelengths (e.g., Nordh et al. 
1982, Harvey et al. 1984, Hurt & Barsony 1996) and also 
at longer wavelengths (Casali et al. 1993, Testi & Sargent 
1998, Davis et al. 1999). In particular at higher spatial 
resolution, this emission region breakes up into numerous 
point-like sources. Because of this source crowding, severe 
source confusion within the Lws-beam can be expected 
at Fir wavelengths. In the next section, we shall attempt 
to disentangle the various contributions from the different 
sources. This is possible, because spatially filling-in and 
overlapping Lws data are in existance (see Tabic |^ and 
Fig.0). 



4.2.1. SEDs of confused sources: Spectro-spatial 
deconvolution of the Lws data 

To search for the Fir sources and to deduce their associ- 
ated SEDs we developed a program for a Maximum Likeli- 



200 



100 



i 



-100 









Fig. 10. The positions observed with the Lws and those of 
the sources included in the Maximum Likelihood Method 
deconvolution. Solid circles correspond to the Lws beams 
and the dashed circle especially marks an Lws observa- 
tion, which is taken at a different time than the grid points 
in the map (see Sect 4.2.3). The filled dots refer to the 
FiR/sub-mm sources for which the deconvolution compu- 
tations were performed - in the wavelength range 50 to 
200 /im - whereas the small open rings identify sources 
which were omitted in the final Fir point source mod- 
elling (see the text) 

hood Method (e.g., Lucy 1974) and applied it to the entire 
set of Lws- map data (Table |l|). 

The first of two basic assumptions is that all sources 
are point-like, which is probably reasonable considering 
the large beam size of the Lws (Hpbw ~ 70") . In sup- 
port of this, SMM 1 was found, in the previous sections. 
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Fig. 11. Upper panels: The Lws spectra of SMM 2, SMM 3 and SMM 4, which are the results of the spectro-spatial 
deconvolution technique described in the text. The solid lines illustrate the blackbody fits to the data which are 
shown as histograms. The error-bars are given by (cf. Appendix A): AFs = J2^ ^Fm / where AFm — AFo-f | 
Fo — ^^Fs\ and where AFq are the errors displayed in Fig.^. Lower panels: Modified blackbody fits over the full 
SED range (10 ~ 10"^ /im). Literature and/or Iso data are coded as in Fig.|| 



to be a point source. Secondly, the exact positions of all 
sources are assumed known and kept fixed. These were 
taken from ISOCAM observations in the mid-infrared (Kaas 
1999) and/or from observations in the sub- millimetre by 
Davis et al. (1999). We then pick an arbitary grid point 
within the field of the Lws-map and try to estimate, at 
that location and for a fixed frequency, the radiative flux 
which originated from all sources in the field. This is done 
at all locations and for all Lws wavelengths in an iterative 
scheme as described in more detail in Appendix A. 

For the restauration of the source spectra, the band- 
width corresponded always to the spectral resolution Rx = 
20. For the wavelength dependent beam pattern of the 
Lws, rim.s, the azimuthal averages of the presently best 
known values were used (see Fig.|). This was done for 
all ten detectors of the Lws and out to ~ 150" from the 
beam centre, which would correspond to about the 5% 



level of a perfectly Gaussian beam. The observational ba- 
sis for this was raster mapping of the point-like planet 
Mars (B.M. Swinyard and C. Lloyd 1999, private commu- 
nication). 

In Fig. |l^, the sources which were included in the mod- 
elling are identified and shown together with the pointings 
of the Lws. These sources are recognized as the strongest 
ones at sub-millimetre (SMM 1 — 11, except SMM 7 and 
10) and/or at mid-infrared wavelengths (SVS 2, SVS 4 and 
SVS 20). The spectra were corrected for background emis- 
sion, which was defined along the northern and western 
edges of the map (see Figs.|| and ||). For these prominent 
sources, convergence was typically achieved after five it- 
erations. In general, the computations were halted after 
visual inspection of the results, but tests were run up to 
1000 iterations demonstrating that the solutions were sta- 
ble. In Appendix B, a demonstration of the algorithm's re- 
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Fig. 12. Same as Fig.|ll| but for SVS 2, SVS 4 and SMM6/SVS20. SVS 2 is actually better modelled as extended 
emission in the far-infrared. This is shown in the panel below by the larger upper limits in the submillimetre, which 
refer to a source filling the Lws beam. For SVS 4 and SVS 20, the dotted lines in the lower panels depict the deduced 
3 a upper limits on the flux XFx in the Lws regime 



liability is provided using known source spectra. It is also 
gratifying to note that the procedure left, as could be ex- 
pected, the spectrum of SMM 1 essentially unaltered, i.e. 
the Maximum Likelihood Method solution for this source 
is entirely within the quoted errors. Yet another check was 
provided by an Lws Guaranteed Time observation half a 
beamwidth south (labelled 'Flow' in Table||), which was 
correctly recovered by our method. Since these observa- 
tions were apart in time by 1 year, this result also indicates 
that no significant FiR-intensity variations have occurred 
on that timescale. 

To properly appreciate the results to be presented be- 
low and summarised in Tables ^ and ||, it might be useful 
to remember that the method attempts to retrieve infor- 
mation which is hidden inside the spatial resolution of the 
observations. Therefore, these solutions represent, at best, 
very likely possibilities of what these source spectra might 
look like. Note that the 30 spectral bins were treated as 



independent observations, so that the fact that the result- 
ing SEDs are continuous and have a meaningful appear- 
ance, lends confidence in the method. There is no a priori 
reason that this should be the case (cf. Eqs.Al to A9). 
However, confirmation must ultimately come from direct 
observation. With this proviso in mind, the results of these 
deconvolution computations are examined for individual 
sources in subsequent sections. 



4.2.2. Spectro-spatial deconvolution results: the 
south-eastern sources 

In the south-east part of the map, the level of confusion 
is high, with many sources contributing to the observed 
emission. A few sources are however dominant, viz. in the 
near infrared SMM 6/SVS 20 and at submillimetre wave- 
lengths SMM 3 and SMM 4, respectively. 
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Fig. 13. Similar as Fig, pi] but for SMM 8, HB 1 and HB2, except that only the upper panels have been included 



SMM 2, 3 and 4: Satisfactory solutions were obtained 
for the sources SMM 3 and SMM 4, whereas the broad 
spectrum of SMM 2 cannot plausibly be represented by 
a modified single-temperature blackbody (see Fig.|ll|). 
This is consistent with the interferometric observations 
by Hogerheijde et al. (1999) and by Choi et al. (1999), 
which failed to reveal a central condensation for SMM 2. 
Kaas (1999) found a K-h&nd nebulosity at this position. 
We interpret its Lws spectrum therefore to be that of an 
extended structure. 

Dust teinperatures of SMM 2 - 4 are of the order of 
30 K (Table^), which is in general agreement with earlier 
work (Harvey et al. 1984), but which again contrasts to the 
lower values obtained by Hurt & Barsony (1996). Also the 
extent of the FiR photospheres is considerably smaller, on 
the order of a few arcsec rather than ten arcsec. For these 
optically thick envelopes (Tab. ||) the masses, which are in 
the range 2 to 6 Mq , would be severly underestimated by 
the optically thin approximations (Tab.^. 

SVS2, SVS4 and SMM6/SVS20: The point source 
assumption leads immediately to inconsistencies for 
SVS 2, which would become unreasonably bright in the far 
infrared, which is highly unlikely given the non-detection 
of the source at longer wavelengths (Fig.|l^). An overall 
better solution is obtained for an extended source, filling 
the Lws beam. The modelling by Hurt & Barsony (1996) 
gave a very different picture, which yielded three nearly 
equally strong Fir point sources, with SVS 2 being only 
slightly stronger at 100 //m. 

The modelling results in largely undetectable Lws 
fluxes for SVS 4 and SVS 20. In particular for SVS 4, these 
low flux levels are consistent with the lack of detection 
at 3 mm, whereas for SVS 20 some emission is possibly 
present at the shortest Lws wavelengths. The objects of 
this group are clearly visible at near infrared wavelengths, 
as stellar clusters and nebulosity and in Table |[ these 
objects are therefore identified as NiR sources. Further, 



our solutions for the Lws data provide reasonable com- 
plements to the CvF data (see Fig.|2|), which have been 
obtained at the considerably higher spatial resolution of 
6" and, hence, are as such not dependent on any models. 

SMM 8, SMM 11 /HB 1 and HB 2: These sources repre- 
sent the fainter members of the south-eastern group in 
the Serpens cloud core. The objects designated PS (Point 
Source) by Hurt & Barsony (1996) are labelled HB 1 and 
HB 2 in Figs. |l0| and 13, Whereas the algorithm finds rea- 
sonable solutions for SMM 8 and HB2, identifying them 
as point sources of roughly 5" size at T>20K, it fails 
for SMMll/HBl. At that location, Kaas (1999) recently 
found a source at 2 /um. 

4.2.3. Spectro-spatial deconvolution results: the 
north-western sources 

The members of this group, north of SMMl, include 
SMM9/S68N, SMM 5 and DEOS, the 'Deeply Embed- 
ded Outburst Star' of Hodapp et al. (1996). As previously 
pointed out, the strong Fir source SMM 1 is comfortably 
displaced from the other three sources, the Lws data of 



which are quite severely confused (see Fig. 10). These con 



fusion problems have also been encountered by Hurt & 
Barsony (1996), who divided the total flux into three equal 
parts among these objects. 

SMM 5 and SMM9/S68N: At the position of SMM 5, 
no recognisable point source was found at Lws wave- 
lengths (Fig.jlJ). The object coincides with a bow-shock 
shaped nebula and it could be that no well defined pre- 
stellar condensation exists. On the other hand, it is quite 
likely that SMM 5 simply is relatively too weak to be de- 
tectable. In order not to unnecessarily increase the noise, 
the source was excluded from further modelling. 

Towards the S 68N complex, two independent Lws ob- 
servations exist: one is a pointed observation, the other a 
grid point in our map. Spatially, these two observations 



B. Larsson et al.: Iso-Lws map of the Serpens cloud core 
SMM 5 SMM 9/ S68N 



13 



DEOS 




80 120 160 40 

Wavelength [/u,m] 



80 120 160 

Wavelength [/-cm] 



80 120 160 

Wavelength [fJ-m] 




10' 10' 10' 

Wavelength [fj-m 



10' 10' 10' 

Wavelength [/U.m] 



10' 10' 10' 

Wavelength [yixm] 



10- 



Fig. 14. Same as Fig.|lT| but for SMM 5, SMM9/S68N and DEOS. The data shown by the asterisks in the SED of 
SMM9/S68N are from Wolf-Chase et al. (1998) 



are separated by about half a beamwidth (33"), and tem- 
porally, by about half a year (Table ^ . The difference in 
flux amounts to about a factor of two, as would be ex- 
pected, if S 68N solely would dominate the Fir emission 
in this region. The observations by Davis et al. (1999) re- 
vealed SMM9/S68N to be the dominant submm source, 
by factors of about 2 - 3 in flux. We expect, therefore, 
this source to be the brightest also in the far infrared, i.e. 
at the longest Lws wavelengths. This seems indeed to be 
the shown in Fig. |lj, where the modelling reveals 

a relatively cool source (T > 20 K). 

DEOS: A further underlying assumption of our decon- 
volution method is that the emission does not vary with 
time. This was previously taken tacitly for granted, since 
this condition is normally fulfilled in the far infrared. For 
the outburst source DEOS, there is a potential risk, how- 
ever, that this may not apply. Between August 1994 and 
July 1995, this object had increased its brightness by al- 
most 5 mag in the X-band (Hodapp et al. 1996). 



At near to midinfrared wavelengths, DEOS was ob- 
served, during about one year and a half, by T. Prusti 
with the Sws on six different dates (see Table|^: April 
14, 1996, to October 22, 1997). The time evolution of the 
short-wave SED is displayed in Fig.|l^, from which it is 
evident that, between April and October 1996, the Sws 
flux had dropped by a factor of two. This is remarkably 
similar to what happened to the (quasi-)simultaneous Lws 
spectra and led us to suspect that some or most of the Fir 
flux in fact was due to DEOS rather than to SMM 9, as 
might have been concluded on the basis of the discussion 
in the previous paragraph. 

Fortunately, the possibility exists to check the tem- 
poral behaviour of the Sws data with independent Iso- 
CAM observations (see Table We performed aperture 
photometry on CAM images, which were obtained in 
two broad-band filters centred on ~7 (LW2) and 15/Ltm 
(LW3), respectively. These observations were performed 
when DEOS was both in the Sws-high (April 14, 1996) 
and in the Sws-low state (September 22, 1997). Table| 
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Fig. 15. Time evolution of the shortwave SED of DEOS, 
obtained with the Sws in the wavelength interval 2 - 
45 ^m. After April 14, 1996, the spectral flux density had 
decreased by a factor of about two and then stayed con- 
stant until the last measurement, i.e. October 22, 1997 

Table 4. Aperture photometry in ISOCAM images in 
Jy/beam, where the beam is that of the Sws 



Object 


Filter 


14 Apr 96 


22 Sep 97 


Ratio 


DEOS 


LW2 


2.8 


2.2 


1.27 


DEOS 


LW3 


6.8 


5.4 


1.26 


EC 53 


LW2 


0.90 


0.90 


1.00 


EC 53 


LW3 


1.12 


1.15 


0.97 



summarises the photometry results for DEOS and a refer- 
ence star (EC 53), expressed in Jy/beam, where the 'per 
beam' refers to the aperture size of the Sws. Evidently, 
the Cam data also reveal significant variability for DEOS, 
by a factor of about 1.3 in both filters. At 7/um, the Sws 
spectra are too noisy to be meaningfully measured, but the 
Sws data in the 15 ^m region are in excellent agreement 
with the CAM-value. 

In conclusion, the evidence presented firmly establishes 
that DEOS had changed its flux in the near- to mid- 
infrared between April and October, 1996. It seems very 
likely therefore that also the Lws spectra show traces of 
this variability. In this context we recall that in the sub- 
millimetre SMM 9/S 68N is the strongest source. Hence, it 
seems very likely that, at the time of the Lws observations 
in October 1996, DEOS was not dominant at the longest 
wavelengths, but might have contributed significantly at 
the shorter Lws wavelengths. A longwave fit for S 68N to 
the pointed Lws observation was properly scaled and sub- 
tracted and the spectro-spatial deconvolution modelling 
was done for the residual. This led finally to the solutions 
for the SEDs of DEOS presented in Fig.|l|. 

The joining of the observed Sws and modelled Lws 
spectra for both the high and the low state is quite spec- 
tacular. This indicates that the algorithm we have used 
is indeed capable of providing satisfactory results. We 
suggest therefore that the spectra presented in Fig. ^ 
represent the most likely SEDs for the outburst source 
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Fig. 16. SEDs of DEOS in high (left) and low state 
(right), referring to the dates shown above the frames. 
The Sws parts of the data have been directly observed, 
whereas the Lws branches result from the spectro-spatial 
deconvolution modelling. Plus-signs refer to data by Ho- 
dapp et al. (1996), and are thus not simultaneous with the 
ISO data, and filled circles to the IsoCAM data presented 
in Table^; other symbols are as in Fig.||. The SEDs con- 
tain considerable spectral structure. Using the calibration 
by Whittet (1998) the optical depth of the silicate feature 
at 10/im, T9.7 ~ 1.3, would implicate some 20 magnitudes 
of visual extinction to the continuum source 
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Fig. 17. Iso observations of FUORs resulted in the in- 
frared SEDs of RNOIB, ZCMa, HH57, V1057Cyg, 
V1331Cyg and V1735Cyg. For reference the SEDs of 
DEOS arc shown to scale 



DEOS. As far as we are aware, these would be seen for 
the first time and should be helpful in constraining out- 
burst models and, perhaps, in differentiating source type 
(e.g. FUOR vs non-FUOR or EXOR). The far-infrared 
SEDs of FUORs appear very different though (Fig. |l^) . 
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4-3. Summarising discussion 

Tabs. I and H summarise the results of the Lws obser- 
vations and our spectro-spatial deconvolution modelling. 
For cool sources, the good wavelength coverage and reso- 
lution of the Lws is generally better adapted to estimate 
the temperature than the poor resolution and short range 
of Iras data. The source emission is generally optically 
thick in the Lws spectral range. Compared to the esti- 
mates based on iRAS data by Hurt & Barsony (1996) our 
temperatures of the SMM sources are generally higher, in- 
dicative of a larger number of photons heating the dust. 
This is also effectively expressed by our larger luminosi- 
ties, although the derived sizes of the FiR photospheres 
(e.g., Tioo^im ~ 1) are generally smaller than the dimen- 
sions obtained by Hurt & Barsony (1996). 

We have deduced the observed f3 values from the slopes 
of the Rayleigh- Jeans part of the source SEDs and obtain 
generally quite shallow distributions, viz. |/3|< 1. This is 
in line with the findings by, e.g., Casali et al. (1993), Testi 
& Sargent (1998) and Davis et al. (1999). However, for a 
given source, the scatter among different observations is 
large and it appears doubtful that (3 could be a useful di- 
agnostic for SMM-type objects. In fact, the observed flat- 
ness probably indicates that the assumption of optically 
thin emission is not valid. This hypothesis is supported by 
detailed radiative transfer models, in which the temper- 
ature distribution is self-consistently calculated. Models 
were calculated for five sources, for which sufficient data 
were available, viz. SMM 1, 2, 3, 4 and 9. These models are 
optically thick out to millimetre wavelengths and envelope 
masses are in the range 2-6 Mq . 

The luminosities of these SMMs are in the range 10- 
70 Lq. Much or most of this luminosity is presumably pro- 
duced by mass accretion processes. These depend on the 
centrally accumulated mass. It is at present difficult to tell 
whether the SMMs will end up on the main sequence as 
low-mass or as intermediate mass stars (;5 1-3 M©; the dis- 
tance too comes in here as a crucial parameter). As OH 
and/or H2O masers normally are found towards young 
sources of higher mass, the detection of such emission to- 
wards SMMl by Rodriguez et al. (1989) and Curiel et 
al. (1993), respectively, would speak in favour of the in- 
termediate mass option, at least for that object. In any 
case, the existence of these massive envelopes around the 
SMMs argues for the idea that these sources are in a very 
early stage of their development. One can also speculate 
that most of these massive envelopes will have become 
dispersed on a relatively short timescale (~ 10^ yr). 

5. Conclusions 

The main conclusions of this work can be summarised as 
follows: 

• In order to disentangle the contribution of individual 
sources to the Fir emission in a 8' x 8' map obtained 



with the Iso-Lws a Maximum Likelihood Method is 
introduced. This spectro-spatial deconvolution tech- 
nique enabled us to restore the Fir spectra (50 to 
200 fim) of previously known submillimetre sources 
(SMMs). 

• Observations of SMMs with the Lws are advantagous 
as the data sample the peak of the spectral energy dis- 
tributions (SEDs). This permits an accurate determi- 
nation of the dust temperatures. These temperatures 
are generally found to be higher than those obtained 
from high resolution Iras data. 

• Fits to the observed SEDs by modified blackbodies 
reveal that the SMMs are generally optically thick at 
Lws wavelengths. 

• In addition, for reasonable assumptions about the 
grain opacities, the Rayleigh- Jeans part of the ob- 
served SEDs is significantly fiatter than what would 
be expected for optically thin dust emission. We inter- 
pret this to indicate the SMMs to be optically thick 
out too much longer wavelengths than previously as- 
sumed. 

• Self-consistent radiative transfer calculations confirm 
the correctness of this assertion. These models of the 
SMMs are optically thick out to mm-wavelengths. 
These dust envelopes are massive, several Mq , suggest- 
ing that these sources are still in their very infancy. 

• The outburst source DEOS has been observed with 
various instruments aboard ISO and at different times. 
The source had certainly varied in infrared brightness 
and we present the full infrared spectra, at both high 
and low states, for the first time. 
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Appendix A: Spectro-spatial deconvolution: the 
algorithm 

The aim of the presented algorithm is to restore the SEDs of 
sources in a crowded field. To be able to handle spatially under- 
sampled and noisy spectroscopic data we developed a program 
based on the Maximum Likelihood Method (e.g., Lucy 1974). 
The present discussion is based on our Lws observations, but 
the method per se should of course be more widely applicable. 

The first of two basic assumptions is that all sources are 
point-like, which is probably reasonable considering the large 
beam size of the Lws (Hpbw~ 70"). Secondly, the exact posi- 
tions of all sources are assumed known and kept fixed. We then 
pick an arbitary grid point within the field of the Lws- map and 
try to estimate, at that location and for a fixed frequency, the 
radiative fiux which originated from all sources in the field. 
This is done at all locations and for all Lws wavelengths in an 
iterative scheme as described in the following. 
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The photon distribution can be assumed to follow Poisson 
statistics and, hence, the probability to observe Om detector 
counts at the position m can be written as 



p — 



(A.l) 



where 

Ns 

E m — / ^m,s Ss 



(A.2) 



is the expected number of counts at position m, arriving 
from Ns sources, each contributing Ss counts and being at 
positions s. Q.m,s is the diffraction pattern of the instrument 
(the 'beam', see Fig.^. Designating the likelihood to observe 
the counts from a map containing Nm points by 



L^PiPiPa ■■■ Pnm 
one can also write 

In L = ^ Om In Em - - In Om! 

m — 1 

Solving for the maximum likelihood, viz. 
OlnL ^ 



(A.3) 



(A.4) 



dSs 



(A.5) 



for all positions s, or 



(A.6) 



From the last expression, a correction factor can be ob- 
tained 



m=l F^/ 
Z^m = l 



(A.7) 



where Fm is now the observed flux (for a linear detector) 
at position m and where 



Ns 

pM _ V n 



(A.8) 



is the expected flux at this position, originating from Ns 
sources at positions s and contributing the flux Ff . Finally, 
the iteration is then 



F^ {New) = C,F/ (Old) 



for all wavelengths A, yielding the source spectra. 



(A.9) 
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Fig. B.l. The Fir spectra of the four point sources in- 
cluded in the test example have simple forms in the wave- 
length range of the Lws observations. The binning corre- 
sponds, as before, to R\ = 20 



200 



100 



i 



-100 



O 



I I I I I I I I I I I I I I I I 







(J 





0OO 





I I 



200 



100 

AR.A. [arcsec] 



-100 



-200 



Fig. B.2. Map showing the grid points of the simulated 
observation, the Hpbw of the Lws and the positions of the 
point sources 1 to 4 having the spectra shown in Fig. DB.ll 



Appendix B: Spectro-spatial deconvolution: test 

As an example of a test case, an observation of the Lws map 
was simulated. We chose four Fir spectra, the shape of which is 
described by simple geometric flgures (positi ve ra mp, negative 
ramp, rectangular box and triangle: see Fig. B.l). The point- 



like sources of these spectra wer e dis tributed on the Lws grid 
in the manner depicted in Fig. [B.2| , i.e. source 1 centred on 
a grid point (which therefore should not be affected by the 
convolution, but still by the noise) , source 2 on the half-power 
contour of the Lws beam, source 3 situated half-way between 
two grid points and source 4 farthest away from any grid point. 

At each wavelength (spectral resolution R\ = 20) for every 
map position the sources were convolved with the Lws beam 
of Fig. ^ and (white) noise was added. This results in the sim- 
ulated observations shown in Fig. B.2 . These "observed" maps 
(one per spectral bin) were then run thrugh the Maximum 
Likelihood algorithm for a number of iterations. 



The result of the deconvolution is shown in Fig. B.4, where 
the succession of iterations is ordered vertically for each of the 
four sources. The flrst row shows the start conditions, where 
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Fig. B.3. The simulated observation of sources 1 to 4 
(with added noise) resuhs in the shown undersampled Lws 
map of spectra (cf. Fig.||). Notice in particular the very 
faint and noisy appearance of source 4 



Fig. B.4. The iterative spectro-spatial deconvolution pro- 
cess is shown vertically for the sources 1 to 4, from left 
to right. The first row shows the start conditions for the 
source spectra, viz. simply a constant spectral distribu- 
tion. Convergence is achieved already after four iterations, 
shown in the fourth row. Iterating further does not im- 
prove on the results: see rows 5 and 6, showing the results 
after 8 and 16 iterations, respectively 



we "guessed" a simple constant spectrum in all cases. These 
straight line spectra have changed to resemble the true source 
spectra already after the first run. In this particular test case, 
the deconvolutions have converged after about 4 iterations 
(fourth row), i.e. successive iterations do not give further im- 



provements. Apparently, the spectra of particularly source 1, 
as expected, but also of sources 2 and 3 are very well r epro - 
duced. The most difficult test case is source 4 (cf also Fig. B.3), 
which is located farthest away from the map points, and there- 
fore considerably more affected by the noise. Nevertheless, our 
method reproduced its spectrum reasonably well too. 

We conclude that this test procedure was successful, which 
lends credibility to the reliability and potential of the devel- 
oped spectro-spatial deconvolution algorithm. 
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